DEM is the English abbreviation of Digital Elevation Model, which is the important original data of watershed topography and feature recognition.DEM is based on the principle that the watershed is divided into cells of m rows and n columns, the average elevation of each quadrilateral is calculated, and then the elevation is stored in a two-dimensional matrix.Since DEM data can reflect local topographic features with a certain resolution, a large amount of surface morphology information can be extracted through DEM, which includes slope, slope direction and relationship between cells of watershed grid cells, etc..At the same time, the surface flow path, river network and watershed boundary can be determined according to certain algorithm.Therefore, to extract watershed features from DEM, a good watershed structure pattern is the premise and key of the design algorithm. Elevation data map 1km data formed according to 1:250,000 contour lines and elevation points in China, including DEM, hillshade, Slope and Aspect maps. Data set projection: Two projection methods: Equal Area projection Albers Conical Equal Area (105, 25, 47) Geodetic coordinates WGS84 coordinate system
TANG Guoan
This data set comprises the plateau soil moisture and soil temperature observational data based on the Tibetan Plateau, and it is used to quantify the uncertainty of model products of coarse-resolution satellites, soil moisture and soil temperature. The observation data of soil temperature and moisture on the Tibetan Plateau (Tibet-Obs) are from in situ reference networks at four regional scales, which are the Nagqu network of cold and semiarid climate, the Maqu network of cold and humid climate, and the Ali network of cold and arid climate,and Pali network. These networks provided representative coverage of different climates and surface hydrometeorological conditions on the Tibetan Plateau. - Temporal resolution: 1hour - Spatial resolution: point measurement - Measurement accuracy: soil moisture, 0.00001; soil temperature, 0.1 °C; data set size: soil moisture and temperature measurements at nominal depths of 5, 10, 20, 40 - Unit: soil moisture, cm ^ 3 cm ^ -3; soil temperature, °C
BOB Su, YANG Kun
The past frozen soil map of the Tibetan Plateau was based on a small number of temperature station observations and used a classification system based on continuity. This data set used the geographically weighted regression model (GWR) to synthesize MODIS surface temperature, leaf area index, snow cover ratio and multimodel soil moisture forecast products of the National Meteorological Information Center through spatiotemporal reconstruction. In addition, precipitation observations of more than 40 meteorological stations, the precipitation products of FY2 satellite observations and the multiyear average temperature observation data of 152 meteorological stations from 2000 to 2010 were integrated to simulate the average temperature data of the Tibetan Plateau, and the permafrost thermal condition classification system was used to classify permafrost into several types: Very cold, Cold, Cool, Warm, Very warm, and Likely thawing. The map shows that, after deducting lakes and glaciers, the total area of permafrost on the Tibetan Plateau is approximately 1,071,900 square kilometers. Verification shows that this map has higher accuracy. It can provide support for future planning and design of frozen soil projects and environmental management.
RAN Youhua, LI Xin
The China Meteorological Forcing Dataset (CMFD) is a high spatial-temporal resolution gridded near-surface meteorological dataset that was developed specifically for studies of land surface processes in China. The dataset was made through fusion of remote sensing products, reanalysis dataset and in-situ observation data at weather stations. Its record starts from January 1979 and keeps extending (currently up to December 2018) with a temporal resolution of three hours and a spatial resolution of 0.1°. Seven near-surface meteorological elements are provided in CMFD, including 2-meter air temperature, surface pressure, specific humidity, 10-meter wind speed, downward shortwave radiation, downward longwave radiation and precipitation rate.
YANG Kun, HE Jie, WENJUN TANG , LU Hui, QIN Jun , CHEN Yingying, LI Xin
Qinghai Tibet Plateau is the largest permafrost area in the world. At present, some permafrost distribution maps have been compiled. However, due to the limited data sources, unclear standards, insufficient verification and lack of high-quality spatial data sets, there is great uncertainty in drawing Permafrost Distribution Maps on TP. Based on the improved medium resolution imaging spectrometer (MODIS) surface temperature (LSTS) model of 1 km clear sky mod11a2 (Terra MODIS) and myd11a2 (Aqua MODIS) product (reprocessing version 5) in 2003-2012, the data set simulates the distribution of permafrost and generates the permafrost map of Qinghai Tibet Plateau. The map was verified by field observation, soil moisture content and bulk density. Permafrost attributes mainly include: seasonally frozen ground, permafrost and unfrozen ground. The data set provides more detailed data of Permafrost Distribution and basic data for the study of permafrost in the Qinghai Tibet Plateau.
ZHAO Lin
The dataset includes soil physical and chemical attributes: pH value, organic matter fraction, cation exchange capacity, root abundance, total nitrogen (N), total phosphorus (P), total potassium (K), alkali-hydrolysable N, available P, available K, exchangeable H+, Al3+, Ca2+, Mg2+, K+ , Na+, horizon thickness, soil profile depth, sand, silt and clay fractions, rock fragment, bulk density, porosity, structure, consistency and soil color. Quality control information (QC) was provided. The resolution is 30 arc-seconds (about 1 km at the equator). The vertical variation of soil property was captured by eight layers to the depth of 2.3 m (i.e. 0- 0.045, 0.045- 0.091, 0.091- 0.166, 0.166- 0.289, 0.289- 0.493, 0.493- 0.829, 0.829- 1.383 and 1.383- 2.296 m) for convenience of use in the Common Land Model and the Community Land Model (CLM). 1.THSCH.nc: Saturated water content of FCH 2.PSI_S.nc: Saturated capillary potential of FCH 3.LAMBDA.nc: Pore size distribution index of FCH 4.K_SCH.nc: Saturate hydraulic conductivity of FCH 5.THR.nc: Residual moisture content of FGM 6.THSGM.nc: Saturated water content of FGM 7.ALPHA.nc: The inverse of the air-entry value of FGM 8.N.nc: The shape parameter of FGM 9.L.nc: The pore-connectivity parameter of FGM 10.K_SVG.nc: Saturated hydraulic conductivity of FGM 11.TH33.nc: Water content at -33 kPa of suction pressure, or field capacity 12.TH1500.nc: Water content at -1500 kPa of suction pressure, or permanent wilting point
DAI Yongjiu, SHANGGUAN Wei
China's second glacier inventory uses the high-resolution Landsat TM/ETM+ remote sensing satellite data as the main glacier boundary data source and extracts the data source with the latest global digital elevation model, SRTM V4, as the glacier attribute, using the current international ratio threshold segmentation method to extract the glacier boundary in bare ice areas. The ice ridge extraction algorithm is developed to extract the glacier ice ridge, and it is used for the segmentation of a single glacier. At the same time, the international general algorithm is used to calculate the glacier attributes, so that the vector data and attribute data that contain the glacier information of the main glacier regions in west China are obtained. Compared with some field GPS field measurement data and higher resolution remote sensing images (such as from QuickBird and WorldView), the glacial vector data in the second glacier inventory data set of China have higher positioning accuracy and can meet the requirements for glacial data in national land, water conservancy, transportation, environment and other fields. Glacier inventory attributes: Glc_Name, Drng_Code, FCGI_ID, GLIMS_ID, Mtn_Name, Pref_Name, Glc_Long, Glc_Lati, Glc_Area, Abs_Accu, Rel_Accu, Deb_Area, Deb_A_Accu, Deb_R_Accu, Glc_Vol_A, Glc_Vol_B, Max_Elev, Min_Elev, Mean_Elev, MA_Elev, Mean_Slp, Mean_Asp, Prm_Image, Aux_Image, Rep_Date, Elev_Src, Elev_Date, Compiler, Verifier. For a detailed data description, please refer to the second glacier inventory data description.
LIU Shiyin, GUO Wanqin, XU Junli
This data set is based on China's second inventory data, Landsat series optical image data with a spatial resolution of 30 meters and cloud coverage of less than 10% and SRTM and other data using ArcGIS, ENVI, Google Earth and other processing software and extracting the glacial lake boundary within 10 km of the glacier boundary by artificial visual interpretation. In addition, the data set adds attributes such as glacial lake type, the mountain range, the province, and the basin to the data as well as quality checking and accuracy verification for the interpreted data. The spatial resolution is 30 meters. It consists of two parts: the glacial lake distribution area vector file and the Inventory Data set of glacial lakes in west China in 2015. It can provide reference data for glacial lake-glacier coupling, water resource utilization and management in west China and can also be used as basic data for regional climate change and cryospheric studies.
WANG Xin
This dataset uses daily temperature data from SMMR (1978-1987), SSM/I (1987-2009) and SSMIS (2009-2015). It is generated by the dual-index (TB, 37v, SG) freeze-thaw discrimination algorithm. The classification results include the frozen surface, the thawed surface, the deserts and water bodies. The data coverage is the main part of China’s mainland, with a spatial resolution of 25.067525 km via the EASE-Grid projection method, and it is stored in ASCIIGRID format. All the ASCII files in this data set can be opened directly with a text program such as Notepad. Except for the head file, the body content is numerically characterized by the freeze/thaw status of the surface soil: 1 for frozen, 2 for thawed, 3 for desert, and 4 for precipitation. If you want to use the icon for display, we recommend using the ArcView + 3D or Spatial Analyst extension module for reading; in the process of reading, a grid format file will be generated, and the displayed grid file is the graphical expression of the ASCII file. The read method comprises the following. [1] Add the 3D or Spatial Analyst extension module to the ArcView software and then create a new View. [2] Activate View, click File menu, and select the Import Data Source option. When the Import Data Source selection box pops up, select ASCII Raster in the Select import file type box. When the dialog box for selecting the source ASCII file automatically pops up, click to find any ASCII file in the data set, and then press OK. [3] Type the name of the Grid file in the Output Grid dialog box (it is recommended that a meaningful file name is used for later viewing) and click the path to store the Grid file, press OK again, and then press Yes (to select integer data) and Yes (to put the generated grid file into the current view). The generated files can be edited according to the Grid file standard. This completes the process of displaying an ASCII file into a Grid file. [4] In the batch processing, the ASCIGRID command of ARCINFO can be used to write AML files, and then use the Run command to complete the process in the Grid module: Usage: ASCIIGRID <in_ascii_file> <out_grid> {INT | FLOAT}. The production of this data is supported by the following Natural Science Foundation Projects: Environmental and Ecological Science Data Center of West China (90502010), Land Data Assimilation System of West China (90202014) and Active and Passive Microwave Radiation Transmission Simulation and Radiation Scattering Characteristics of the Frozen Soil (41071226).
LI Xin
A multi-layer soil particle-size distribution dataset (sand, silt and clay content), based on USDA (United States Department of Agriculture) standard for regional land and climate modelling in China. was developed The 1:1,000,000 scale soil map of China and 8595 soil profiles from the Second National Soil Survey served as the starting point for this work. We reclassified the inconsistent soil profiles into the proper soil type of the map as much as possible because the soil classification names of the map units and profiles were not quite the same. The sand, silt and clay maps were derived using the polygon linkage method, which linked soil profiles and map polygons considering the distance between them, the sample sizes of the profiles, and soil classification information. For comparison, a soil type linkage was also generated by linking the map units and soil profiles with the same soil type. The quality of the derived soil fractions was reliable. Overall, the map polygon linkage offered better results than the soil type linkage or the Harmonized World Soil Database. The dataset, with a 1-km resolution, can be applied to land and climate modelling at a regional scale. Data characteristics: projection:projection Coverage: China Resolution: 0.00833 (about 1 km) Data format: FLT, TIFF Value range: 0%-100% Document describing: Floating point raster files include: Sand1. FLT, clay1. FLT -- surface (0-30cm) sand, clay content. Sand2. FLT, clay2. FLT -- content of sand and clay in the bottom layer (30-100cm). PSD. HDR -- header file: Ncols - the number of columns Nrows- rows Xllcorner - latitude in the lower left corner Yllcorner - longitude of the lower left corner Cellsize - cellsize NODATA_value - a null value byteorder - LSBFIRST, Least Significant Bit First TIFF raster files include: Sand1. Tif, clay1. Tif - surface (0-30cm) sand, clay content. Sand2. Tif, clay2. Tif - bottom layer (30-100cm) sand, clay content.
SHANGGUAN Wei, DAI Yongjiu
Soil bulk density, porosity, water content, water characteristic curve, saturated hydraulic conductivity, particle analysis, infiltration rate, and sampling point location information in the upper reaches of the Heihe River Basin. 1. The data is for 2014 supplementary sampling for 2012, using the ring knife to take the original soil; 2. The soil bulk density is the dry bulk density of the soil and is measured by the drying method. The original ring-shaped soil sample collected in the field was thermostated at 105 ° C for 24 hours in an oven, and the soil dry weight was divided by the soil volume (100 cubic centimeters) , unit: g/cm 3 . 3. Soil porosity is obtained according to the relationship between soil bulk density and soil porosity; 4. Soil infiltration analysis data set, the data is the field experimental measurement data from 2013 to 2014. 5. The infiltration data is measured by “MINI DISK PORTABLE TENSION INFILTROMETER”, and the approximate saturated hydraulic conductivity under a certain negative pressure is obtained. 6. Soil particle size data was measured at the Grain Granulation Laboratory of the Key Laboratory of the Ministry of Education of Lanzhou University. The measuring instrument is a Malvern laser particle size analyzer MS2000. 7. The saturated hydraulic conductivity is measured according to the enamel hair self-made instrument of Yi Yanli (2009). The Marioot bottle was used to maintain the head during the experiment; at the same time, the Ks measured at the time was converted to the Ks value at 10 °C for analysis and calculation. 8. Soil water content data is measured using ECH2O, including 5 layers of soil water content and soil temperature. 9. The water characteristic curve is measured by the centrifuge method: the undisturbed soil of the ring cutter collected in the field is placed in a centrifuge, and each of the speeds is measured at 0, 310, 980, 1700, 2190, 2770, 3100, 5370, 6930, 8200, 11600. The secondary rotor weight is obtained.
HE Chansheng
Snow cover dataset is produced by snow and cloud identification method based on optical instrument observation data, covering the time from 1989 to 2018 (two periods, from January to April and from October to December) and the region of Qinghai-Tibet Plateau (17°N-41°N, 65°E-106°E) with daily product, which takes equal latitude and longitude projection with 0.01°×0.01° spatial resolution, and characterizes whether the ground under clear sky or transparent thin cloud is covered by snow. The input data sources include AVHRR L1 data of NOAA and MetOp serials of satellites, and L1 data corresponding to AVHRR channels taken from TERRA/MODIS. Decision Tree algorithm (DT) with dynamic thresholds is employed independent of cloud mask and its cloud detection emphasizes on reserving snow, particularly under transparency cirrus. It considers a variety of methods for different situations, such as ice-cloud over the water-cloud, snow in forest and sand, thin snow or melting snow, etc. Besides those, setting dynamic threshold based on land-surface type, DEM and season variation, deleting false snow in low latitude forest covered by heavy aerosol or soot, referring to maximum monthly snowlines and minimum snow surface brightness temperature, and optimizing discrimination program, these techniques all contribute to DT. DT discriminates most snow and cloud under normal circumstances, but underestimates snow on the Qinghai-Tibet Plateau in October. Daily product achieves about 95% average coincidence rate of snow and non-snow identification compared to ground-based snow depth observation in years. The dataset is stored in the standard HDF4 files each having two SDSs of snow cover and quality code with the dimensions of 4100-column and 2400-line. Complete attribute descriptions is written in them.
ZHENG Zhaojun, CHU Duo
This data is originated from the 1:100,000 national basic geographic database, which was open freely for public by the National Basic Geographic Information Center in November 2017. The boundary of the Qinghai-Tibet Plateau was spliced and clipped as a whole, so as to facilitate the study on the Qinghai-Tibet plateau. This data set is the 1:100,000 administrative boundaries of the qinghai-tibet plateau, including National_Tibet_line、 Province_Tibet、City_Tibet、County_Tibet_poly and County_Tibet_line. Administrative boundary layer (County_Tibet_poly) property name and definition: Item Properties Describe Example PAC Administrative division code 513230 NAME The name of the County line name Administrative boundary layer (BOUL) attribute name and definition: Item Properties Describe Example GB classification code 630200 Administrative boundary layer (County_Tibet_line) attribute item meaning: Item Properties Describe Example GB 630200 Provincial boundary GB 640200 Prefectural, municipal and state administrative boundaries GB 650201 county administrative boundaries (determined)
National Basic Geographic Information Center
Based on the existing natural hole data of 15 active layer depth monitoring sites in the Qinghai-Tibet Engineering Corridor, the active layer depth distribution map of the Qinghai-Tibet Engineering Corridor was simulated using the GIPL2.0 frozen soil model. The model required synthesis of a temperature data set of time series. The temperature data were divided into two phases according to the time spans, which were 1980-2009 and 2010-2015. The data of the first phase were from the Chinese meteorological driving data set (http://dam. Itpcas.ac.cn/rs/?q=data#CMFD_0.1), and the data of the second phase was the application of MODIS surface temperature products (MOD11A1/A2 and MYD11A1/A2) with a spatial resolution of 1 km. In addition, the soil type data required by the model came from the China Soil Database (V1.1) and have a resolution of 1 km. At the same time, the topography was also considered. The research area was classified into 88 types based on the measured soil thermophysical parameters and land cover types, and then the simulation was performed. The simulation results were compared with the field measured data. The results showed that they were highly consistent, and the correlation coefficient reached 0.75. In alpine areas, the average depth of the active layer is below 2.0 m. However, in the river valleys, the average depth of the active layer is above 4.0 m. In the high plain area, the depth of the active layer is usually between 3.0 m and 4.0 m.
NIU Fujun, YIN Guoan
As the main parameter in the land surface energy balance, surface temperature indicates the degree of land-atmosphere energy and water transfer and is widely used in research on climatology, hydrology and ecology. In the study of frozen soil, climate is one of the decisive factors for the existence and development of frozen soil. The surface temperature is the main climatic factor affecting the distribution of frozen soil and affects the occurrence, development and distribution of frozen soil. It is the upper boundary condition for modelling frozen soil and is significant to the study of hydrological processes in cold regions. The data set was based on the DEM and observation station data of the Tibetan Plateau Engineering Corridor and analysed the changing trend of surface temperature on the Tibetan Plateau from 2000 to 2014. Using the surface temperature data products MOD11A1/A2 and MYD11A1/A2 of MODIS aboard Terra and Aqua, the surface temperature information under cloud cover was reconstructed based on the spatio-temporal information of the images. The reconstruction information and surface temperature representativeness problems were analysed using information obtained from 8 sites, including the Kunlun Mountains (wetland, grassland), Beiluhe (grassland, meadow), Kaixinling (meadow, grassland), and Tanggula Mountain (meadow, wetland). According to the correlation coefficient (R2), root-mean-square error (RMSE), mean absolute error (MAE) and mean deviation (MBE), the following results were obtained: (1) the reconstruction accuracy of MODIS surface temperature under cloud cover is higher when it is based on spatio-temporal information; (2) the weighted average representation is the best when generalizing four observations of Terra and Aqua. By analysing the reconstruction of MODIS surface temperature information and representativeness problems, the average annual MODIS surface temperature data of the Tibetan Plateau and the engineering corridor from 2000 to 2010 were obtained. According to the data set, the surface temperature from 2000 to 2010 also experienced volatile rising trends from 2000 to 2010, which is basically consistent with the changing trend of the climate change in the permafrost regions of the Tibetan Plateau and the Qinghai-Tibet Engineering Corridor.
NIU Fujun, YIN Guoan
The GAME/Tibet project conducted a short-term pre-intensive observing period (PIOP) at the Amdo station in the summer of 1997. From May to September 1998, five consecutive IOPs were scheduled, with approximately one month per IOP. More than 80 scientific workers from China, Japan and South Korea went to the Tibetan Plateau in batches and carried out arduous and fruitful work. The observation tests and plans were successfully completed. After the completion of the IOP in September, 1998, five automatic weather stations (AWS), one Portable Atmospheric Mosonet (PAM), one boundary layer tower and integrated radiation observatory (Amdo) and nine soil temperature and moisture observation stations have been continuously observed to date and have obtained extremely valuable information for 8 years and 6 months consecutively (starting from June 1997). The experimental area is located in Nagqu, in northern Tibet, and has an area of 150 km × 200 km (Fig. 1), and observation points are also established in D66, Tuotuohe and the Tanggula Mountain Pass (D105) along the Qinghai-Tibet Highway. The following observation stations (sites) are set up on different underlying surfaces including plateau meadows, plateau lakes, and desert steppe. (1) Two multidisciplinary (atmosphere and soil) observation stations, Amdo and NaquFx, have multicomponent radiation observation systems, gradient observation towers, turbulent flux direct measurement systems, soil temperature and moisture gradient observations, radiosonde, ground soil moisture observation networks and multiangle spectrometer observations used as ground truth values for satellite data, etc. (2) There are six automatic weather stations (D66, Tuotuohe, D105, D110, Nagqu and MS3608), each of which has observations of wind, temperature, humidity, pressure, radiation, surface temperature, soil temperature and moisture, precipitation, etc. (3) PAM stations (Portable Automated Meso - net) located approximately 80 km north and south of Nagqu (MS3478 and MS3637) have major projects similar to the two integrated observation stations (Amdo and NaquFx) above and to the wind, temperature and humidity turbulence observations. (4) There are nine soil temperature and moisture observation sites (D66, Tuotuohe, D110, WADD, NODA, Amdo, MS3478, MS3478 and MS3637), each of which has soil temperature measurements of 6 layers and soil moisture measurement of 9 layers. (5) A 3D Doppler Radar Station is located in the south of Nagqu, and there are seven encrypted precipitation gauges in the adjacent (within approximately 100 km) area. The radiation observation system mainly studies the plateau cloud and precipitation system and serves as a ground true value station for the TRMM satellite. The GAME-Tibet project seeks to gain insight into the land-atmosphere interaction on the Tibetan Plateau and its impact on the Asian monsoon system through enhanced observational experiments and long-term monitoring at different spatial scales. After the end of 2000, the GAME/Tibet project joined the “Coordinated Enhanced Observing Period (CEOP)” jointly organized by two international plans, GEWEX (Global Energy and Water Cycle Experiment) and CL IVAR (Climate Change and Forecast). The Asia-Australia Monsoon Project (CAMP) on the Tibetan Plateau of the Global Coordinated Enhanced Observation Program (CEOP) has been started. The data set contains POP data for 1997 and IOP data for 1998. Ⅰ. The POP data of 1997 contain the following. 1. Precipitation Gauge Network (PGN) 2. Radiosonde Observation at Naqu 3. Analysis of Stable Isotope for Water Cycle Studies 4. Doppler radar observation 5. Large-Scale Hydrological Cycle in Tibet (Link to Numaguchi's home page) 6. Portable Automated Mesonet (PAM) [Japanese] 7. Ground Truth Data Collection (GTDC) for Satellite Remote Sensing 8. Tanggula AWS (D105 station in Tibet) 9. Syamboche AWS (GEN/GAME AWS in Nepal) Ⅱ. The IOP data of 1998 contain the following. 1. Anduo (1) PBL Tower, 2) Radiation, 3) Turbulence SMTMS 2. D66 (1) AWS (2) SMTMS (3) GTDC (4) Precipitation 3. Toutouhe (1) AWS (2) SMTMS (3 )GTDC 4. D110 (1) AWS (2) SMTMS (3) GTDC (4) SMTMS 5. MS3608 (1) AWS (2) SMTMS (3) Precipitation 6. D105 (1) Precipitation (2) GTDC 7. MS3478(NPAM) (1) PAM (2) Precipitation 8. MS3637 (1) PAM (2) SMTMS (3) Precipitation 9. NODAA (1) SMTMS (2) Precipitation 10. WADD (1) SMTMS (2) Precipitation (3) Barometricmd 11. AQB (1) Precipitation 12. Dienpa (RS2) (1) Precipitation 13. Zuri (1) Precipitation (2) Barometricmd 14. Juze (1) Precipitation 15. Naqu hydrological station (1) Precipitation 16. MSofNaqu (1) Barometricmd 16. Naquradarsite (1)Radar system (2) Precipitation 17. Syangboche [Nepal] (1) AWS 18. Shiqu-anhe (1) AWS (2) GTDC 19. Seqin-Xiang (1) Barometricmd 20. NODA (1)Barometricmd (2) Precipitation (3) SMTMS 21. NaquHY (1) Barometricmd (2) Precipitation 22. NaquFx(BJ) (1) GTDC(2) PBLmd (3) Precipitation 23. MS3543 (1) Precipitation 24. MNofAmdo (1) Barometricmd 25. Mardi (1) Runoff 26. Gaize (1) AWS (2) GTDC (3) Sonde A CD of the data GAME-Tibet POP/IOP dataset cd (vol. 1) GAME-Tibet POP/IOP dataset cd (vol. 2)
MA Yaoming
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS 0931-4967287 poles@itpcas.ac.cnLinks
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved | No.11010502040845
Tech Support: westdc.cn